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Abstract: We discuss the calculation of charged Higgs boson production in asso- 
ciation with a top quark in the MC@NLO framework for combining NLO matrix 
elements with a parton shower. The process is defined in a model independent man- 
ner for wide applicability, and can be used if the charged Higgs boson mass is either 
greater or less than the mass of the top quark. For the latter mass region, care is 
needed in defining the charged Higgs production mode due to interference with top 
pair production. We give a suitable definition applicable in an NLO (plus parton 
shower) context, and present example results for the LHC 
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1. Introduction 

Even in the absence of any definitive experimental evidence, the Higgs mechanism 
of spontaneous symmetry breaking [1, 2, 3, 4, 5] remains a very promising candidate 
to explain the existence of massive weak gauge bosons. This mechanism is imple- 
mented in its minimal version within the Standard Model (SM), i.e. using a single 
SU{2)l Higgs doublet. Yet, more complex scenarios involving additional scalar fields 
remain possible and could display interesting properties such as, for example, new 
sources of CP violation, and/or embedding in more complex models like the Minimal 
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Supersymmetric Standard Model (MSSM). In this context, existing and future high 
energy colhders will have to determine not only if the elusive Higgs particle exists, 
but also if the observed scalar sector is minimal or not. Any conclusive answer to 
the last question strongly relies upon the possibility of observing a charged Higgs 
boson. Indeed, the discovery of such a particle would clearly imply the presence 
of additional non trivial weak multiplet(s) in the scalar potential responsible for 
spontaneous symmetry breaking. 

If the charged Higgs boson is lighter than the top quark, the most promising 
production mode at hadron colliders (see e.g. [6] for a recent review) is through the 
top quark decay t H^h (or i H~b). If it is larger, as suggested by indirect 
measurements such as the b 57 branching ratio [7], the most promising process 
is the direct production in association with a single top quark, which is the focus 
of the present work. Other production modes for the charged Higgs boson, such 
as s-channel single, pair or associated production processes, are less favorable for 
discovery in most models. 

In order to define as precisely as possible an isolation strategy and (if a charged 
Higgs boson happens to be actually observed) in order to render possible the extrac- 
tion of physical parameters such as its coupling to heavy fermions, accurate predic- 
tions are necessary at the fully exclusive level for this channel. For similar Standard 
Model processes involving a single top quark (which can play the role of backgrounds 
to the considered channel), the current state of the art is the combination of NLO 
parton-level matrix elements with Monte Carlo event generators. Those generators 
use parton shower algorithms to simulate the effect of further soft and coUinearly 
enhanced radiation, as well as modelling hadronization effects and the underlying 
event. One such approach for combining NLO matrix elements with parton showers 
is the MC@NLO algorithm of [8]. All of the three single top production modes have 
all already been implemented in this framework [9, 10], including angular correla- 
tions using the method described in [11] (for a recent study of angular correlations 
in top quark production, see [12]). It is thus natural to implement the production of 
a charged Higgs boson (i/^) with an accompanying single top quark, first calculated 
at NLO in [13, 14]. Single top production in both the s- and i-channel modes was 
also recently implemented in the POWHEG framework for combining NLO matrix 
elements with parton showers [15]. 

As we will see, H'^t production is theoretically similar to the Wt mode (essen- 
tially, one replaces the W boson in all Wt Feynman diagrams with a charged Higgs 
boson). This creates an extra motivation for studying H"^ production in a frame- 
work which combines NLO matrix elements with parton showers. Indeed, it is well 
known that Wt production at NLO mixes with top pair production at LO, followed 
by decay of one of the final state top particles. Thus, the meaning of the Wt mode 
itself, and by analogy the if^ mode when mH± < irit, becomes a matter of careful 
definition, which must be applicable in an experimental context. This then demands 
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Figure 1: Leading order diagrams for single top production (double line) in association 
with a charged Higgs boson (dashed line). 

calculations which are at least as complex as MC@NLO, in that the problem occurs 
only at NLO and beyond, and only in a fully exclusive, all-orders computation can a 
suitable definition of the considered mode be tested within an experimental context. 

The structure of the paper is as follows. In section 2 we outline the computation 
of the charged Higgs production mode at NLO in QCD, together with its imple- 
mentation in the MC@NLO framework. In section 3, we describe how to define the 
H~t mode for niu- < uit in two different ways (whose comparison measures inter- 
ference contributions), leaning on the preceding discussion of Wt production in [10]. 
In section 4 we present example results, also discussing the transition region where 
tuh- — TTit. We conclude in section 5. 

2. H't production at NLO 

In this section we describe our calculation of the H~t process up to NLO in QCD 
perturbation theory. In order to check our results, we performed the calculation using 
two different methods for dealing with infrared divergences in the real and virtual 
contributions. Firstly, we used the Catani-Seymour algorithm of [16, 17]. Given that 
this calculation has not been performed before, we present some salient details in 
what follows. Secondly, we used the FKS subtraction formahsm of [18, 19]. This is 
the formalism required for the implementation of NLO results within MC@NLO, and 
as such has already been used in previous processes. Hence, the details are extremely 
similar to the case of Wt production considered in [10], and we do not present many 
of them here. We start by discussing the LO result, which will also be useful in 
introducing notation, in the following section. 

2.1 Born computation 

The LO Feynman diagrams for charged Higgs production are shown in figure 1. From 
now on, we consider explicitly production (the case of production being the 
same, as shown for the Wt case in [10]). Note that, as in [10], we assume a five fiavor 
scheme [20] for the quark sector, where the b quark is massless. The case where b 
quarks are generated explicitly by initial state gluon splitting was considered in [14]. 
We label momenta as follows: 
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bipi) + gip2)^t{h) + H-{k2). 



(2.1) 



One must define the coupling of the charged Higgs boson to fermions, which is 
clearly a matter of convention. It is of most practical use to leave this coupling 
model- independent, and we choose to write the Higgs-fermion vertex as follows: 

GH-UD^iVuD[AuD-BuDl5], (2.2) 

where U and D are up- and down- type quarks respectively, and we have explicitly 
factored out the corresponding CKM matrix element Vu^,. To simphfy notation from 
now on, we consider the case where the H~ couples only to a i and a b quark (by far 
the most dominant process), and write 

GH-tb = ^Vtb[a-b^,], (2.3) 

where a and b may be complex in general. When presenting results later in the paper, 
we will use the specific example of a type-II two-Higgs doublet model, in which a 
and b are given as in eq. (A. 14) of appendix A, and depend explicitly on the top 
and bottom quark masses. In the matrix element, all quarks other than the top are 
treated as massless, including the b quark. In both the Born and NLO calculations, 
the pure scalar and pure pseudoscalar couplings add incoherently i.e. the total cross- 
section is proportional to|ap-|-|6p. A comment is in order regarding the use of 
the massless b quark approximation. The reader may worry about the fact that nib 
dependence is kept in the Yukawa couplings, but neglected in the matrix elements. 
In actual fact there is no inconsistency, as discussed clearly in [14, 21]. Keeping the 
b mass in the Yukawa coupling merely corresponds to keeping only the leading rrib 
behavior. 

2.2 NLO computation 

As stated above, we performed the calculation at NLO accuracy using two different 
subtraction formalisms for dealing with infrared divergences in order to check our 
results. Before discussing these in more detail, some remarks are in order regarding 
the treatment of ultraviolet divergences in the couplings and heavy quark mass. 

2.2.1 Coupling and mass renormalization 

We evaluated all one-loop diagrams using dimensional regularization in = 4 — 2e 
dimensions. Up to NLO, one finds both double and single poles in e, arising from 
infra-red (IR) and ultraviolet (UV) divergences. The former cancel in the sum of 
virtual and real corrections (or, in the case of initial state coUinear divergences, arc 
removed by counterterms) , as described in the next section. UV divergences are 
removed by renormalization of the strong and Yukawa couplings, and of the top 
quark mass. 
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As in [10], we modify the MS-scheme QCD coupling such that the top quark 
loop contribution is subtracted on-shell. In this scheme [22], the top quark virtual 
contributions decouple in the limit of small external momentum. Specifically, one 
has 
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(2.4) 

where ^ip and ^lr are the factorization and renormalization scales respectively. Fur- 
thermore, /3o = (IIC'a — 2n/)/3, with n/ the number of light flavors (here five) plus 
one. The QCD coupling then satisfies: 
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so that the top quark loop contribution is indeed removed. 

The renormalization of the top quark mass in the on-shell scheme is given by 
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(2.6) 



The renormalization of the Yukawa coupling y (which represents either a or 6 in 
eq. (2.3)) to a quark of mass m is related to the appropriate mass counterterm via 
(see appendix A) 

5y=^, (2.7) 



m 



where the renormalized quark mass appears in the denominator on the right-hand 
side. Note that this result is independent of the mass (since 5m oc m), and also that 
the mass counterterms appearing in the renormalization of both the heavy quark 
mass and Yukawa couplings must be in the same renormalization scheme (in our 
case the on-shell scheme). 



2.2.2 Catani-Seymour subtraction 

Here we briefly summarize the NLO calculation of the H~t process in the Catani- 
Seymour dipole formalism of [16, 17]. Our terminology and notation is similar to 
those papers, to which we refer the reader for more details. 

In the dipole formalism, as in any other subtraction formalism, the real and 
virtual contributions are dealt with so as to render them separately finite and thus 
numerically integrable. This is achieved through a subtraction of the leading phase- 
space divergences in the former, and of the infrared poles in the latter; we shall 
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give more details in the following. The total NLO partonic cross-section has the form 

Jo 



+ / dX2a^''''{x2]PuX2P2]liF). (2-8) 

Jo 

where we have explicitly denoted the dependence on the initial momenta p^, and 
is the factorization scale. The first term on the right-hand side of cq. (2.8) has 2 — > 2 
kinematics and is comprised of the virtual corrections and subtraction term. The 
second term has 2-^3 kinematics, and includes the regularized real corrections. 
The final two terms constitute a finite remainder left after factorization of initial 
state coUinear singularities, and involve integrals over the longitudinal momentum 
fractions Xi of the initial state partons. 

In more detail, the regularized virtual corrections have the form 

<7^^°('Hpi,P2) = j ci$(^HPi,P2)[^2Re [A4i..oopA^Ln 

+ {t,H-g,h\m%H-g,h)\ , (2.9) 

Je=0 

where Ma is the number of color degrees of freedom associated with incoming parton 

a, and 0?$^^-* the phase-space of the two final state particles (including spin averaging). 
Furthermore, the second term denotes the subtraction operator I, which is matrix- 
valued in color space and acts on the color vectors |t, H\ g, h) associated with the 
particles entering the Born interaction. The subtraction operator may be further 
decomposed as 

1(e) = I2(e,//^;mt,/c2) +Ib{e, ii'^;k2,mt,pi) +Ig{e, ii'^;k2,mt,p2) +Ibg{e, ii'^;pi,p2), 

(2.10) 

where !„ includes the contribution from a gluon exchanged (across the final state cut 
in the squared amplitude) between parton a and the top quark, 1;,^ that from a gluon 
exchanged between the initial state partons, and I2 from gluons exchanged between 
the final state particles (this is zero for the process considered here, given that the 
Higgs boson is a color singlet) . 

The real emission term in eq. (2.8) has the schematic form 



jirjir\M3{kuk2,k3;pi,p2)f - ^ T>{ki,k2,ks;pi,p2) 

^ ^ dipoles 



(2.11) 

where is the momentum of the extra parton at NLO, d^^^^ the phase-space of the 
final state particles, and AI3 the real emission amplitude consisting of the following 
four subprocesses: 

(a): b{pi) + g{p2) ^ t{ki) + H-{k2) + g{h)- 
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(b) : g{pi) + g{p2) ^ t{ki) + H-{k2) + 

(c) : q/q{pi) + 6(^2) ^ t{ki) + //"(/cs) + ^/^(/ca); 

(d) : + g(p2) ^ t{k^) + ^-(A;2) + &(A;3). 

Here q denotes a generic quark species (which may be a 6 quark). Note, however, 
that process (a) can be obtained from (b) by crossing. Also, processes (c) and (d) 
interfere when q — b, and care must be taken such that this is correctly dealt with. 
The contribution from these processes is in any case negligible in practice, due to the 
smallness of the b quark parton density. When q ^ b, process (d) is finite as e ^ 0. 

The second term on the right-hand side of eq. (2.11) contains a sum over Catani- 
Seymour dipoles, representing all possible gluon exchanges across the final state cut, 
and collinear splittings of initial state particles. Using the standard nomenclature 
of [16, 17] with upper and lower indices denoting initial and final state particles 
respectively, one may write these as 1)°'^^''^ and T>^^ for an emitter a emitting particle 
P, and 7 the spectator particle. There are also dipoles Vj^ corresponding to final 
state emitters and initial state spectators. For the present process there are no dipoles 
involving final state emitters and spectators, due to the final state color singlet Higgs 
particle. A complete classification of dipoles is then 

. Process (a): Edipoies = + T^'"'' + + + % + 
• Process (b): Edipoies = + + + ^ 

. Process (c): ^^^^^^^^ = P^^-f + , 

where gi^2 denote the two initial state gluons in process (b). 

Finally, the last term of eq. (2.8) (i.e. the finite remainder after collinear factor- 
ization) has the general expression 

[ dxia^^^^^\xi]Xipi,p2]i^l) = [ dxx j d^^^\xipi,p2) 
Jo J 

{ki,k2;xipi,p2 ^(k''''{xi) + P'''' {xipi,x, fip)^ ki,k2;xipi,p2). 

(2.12) 

for the initial state parton whose momentum is pi, with a similar expression for p2. In 
this formula, i labels the parton emerging from the hadron, and i' a particle from the 
collinear splitting which interacts with the second incoming particle. The functions 
K*'*' and P*'*' (matrix- valued in color space and evaluated between color vectors) are 
process- independent, and involve products of ingoing and outgoing momenta. 

2.2.3 FKS subtraction 

The FKS calculation is extremely similar to that carried out for Wt production 
in [10], and thus we do not repeat the discussion here (the calculation is also similar 
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for both rriH- > rrit and mn- < rrit). We used FORM [23] for both the virtual and 
real corrections. 

The virtual corrections in our two calculations were compared analytically, and 
found to agree exactly. The real corrections were compared numerically, and agree- 
ment found within numerical uncertainties. As a further check, we have compared 
the numerical implementation of our total NLO results with that of [14], which is 
available in the pubhcly available Prospino2.0 package [24, 25]. Such checks are pos- 
sible only for rnn- > trit, where one does not have to worry about interference with 
top pair production. For the total NLO cross-section, we find agreement within a 
percent for a range of charged Higgs boson masses between 300 GcV and 1 TeV. 
Note that in order to compare with the calculation of [14], we changed the rcnormal- 
ization of the Yukawa coupling to the MS scheme rather than the on-shell scheme, 
as detailed in that paper. Checks were also performed ior /ip fiji (to exphcitly test 
terms involving logarithms of the form logdip/ fJ-R)), and a similar level of agreement 
found. 

Having checked the implementation of our NLO calculation, the FKS code was 
then interfaced with MC@NLO as prescribed in [8]. Note that there is a techni- 
cal subtlety regarding the collinear limits for partonic subprocesses involving gluon 
branchings. In the case of the FKS formalism, one introduces additional universal 
splitting kernels (as described in [18] and applied also in [10]). See appendix B for 
details. 

2.3 Implementation in MC@NLO 

In this section, we describe how the calculation of the previous section can be im- 
plemented in the MC@NLO framework for combining NLO matrix elements with 
a parton shower algorithm. Given that several processes have already been imple- 
mented in this formalism, we refer the reader to [8, 9, 26] for more technical details, 
and here briefly describe only those aspects that arc relevant for H^t production. 
Furthermore, the implementation of the present process is, as remarked in the pre- 
vious sections, extremely similar to that of Wt production considered in [10] . 

The MC@NLO algorithm first presented in [8] gives a systematic procedure for 
combining a NLO matrix element with a parton shower, in such a way that dou- 
ble counting of radiation is avoided. This is achieved through the definition of the 
so-called Monte Carlo subtraction terms, that are designed to cancel exactly the 
contributions at ncxt-to-leading order to the cross section of interest, given by the 
parton shower. The MC subtraction terms are factorized into process-independent 
kernels, that depend solely on how the Monte Carlo treats the collinear and soft emis- 
sions (e.g. through the definition of the shower variables), times process-dependent 
short-distance cross sections. The latter essentially coincide with Born matrix el- 
ements, and are therefore available as part of the NLO parton-level cross sections. 
As stressed above, the kernels are MC-specific; for a given MC, only a handful of 
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cases have to be considered, corresponding to initial- or final-state branchings, from 
a massless or a massive particle. All relevant computations have been carried out for 
the case where the MC is Fortran HERWIG (see in particular refs. [8, 26, 9]). This is 
now also the case for Herwig++ [27]^; kernels relevant to initial-state radiation have 
also been computed for Pythia [31]. The subtraction terms have the schematic form 

da ^EEE^'^f''^ ' (2-13) 

MC * MC ^ ^ 

i L I 

where i labels different partonic subprocesses, L the leg from which the extra parton 
is emitted (i.e. the parton appearing at NLO which is double counted by the parton 
shower), and / denotes a given color structure. Furthermore, each partonic branching 
has a shower energy scale 

E^, = \qa-qf,\, (2.14) 

where qa and denote the 4-momenta of the color partners, one of which undergoes 
the branching. Thus, the shower scales are completely determined by the partonic 
subprocesses involved in the hard scattering. 

Crucially, the partonic subprocesses and color structures involved in the H~t 
production process are the same as in Wt production. Thus, the subtraction terms 
and shower scales for the present case do not have to be recalculated, and can be 
read off from section 3.3 of [10] (with the W boson replaced by a charged Higgs). 
We refer the reader to that paper for details. 

A given partonic subprocess may have more than one color flow, and it is then 
necessary to select one to feed to the parton shower Monte Carlo. We do this anal- 
ogously to what was used in [10], which can be summarized as follows. Where more 
than one color flow exists, we select one on a statistical basis, with each color flow 
weighted by its large Nc approximation. 

Note that in the present implementation, spin correlations have not yet been 
included in the decay of the top quark. These will be included at a later stage, using 
the technique of [11]. 



3. Interference with tt production 

In the previous subsections, we have discussed the calculation of the H~t production 
process up to next-to-leading order in QCD, and its subsequent implementation in 
the MC@NLO framework for interfacing with a parton shower. In this discussion 
we ignored the fact that for uih- < rrit ei theoretical problem arises in that, at NLO 
and beyond, the single top process interferes with the production of a top quark 
pair, with decay of the antitop particle into a charged Higgs boson and b quark. 

^Some processes interfaced with Herwig-I— I- have already appeared in pubUc MC@NLO codes [28, 
29, 30]. 
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Figure 2: Subset of NLO real emission contributions to H^t production, consisting of top 
pair production with decay of the antitop quark to produce a charged Higgs boson and b 
quark. 



We consider this issue, including how to recover a well-defined meaning of the H t 
process, in this section. Our discussion leans heavily on the results obtained for the 
Wt process in [10, 32]. 

The interference problem can be appreciated by considering figure 2, which shows 
a subset of the NLO real-emission corrections to H~t production. The contribution 
from these diagrams grows as the invariant mass of the final state H~b tends towards 
the top mass rrit, due to the propagator for the intermediate t quark which is moving 
on-shell. In practice this means that the NLO correction to the LO H~t inclusive 
cross-section is huge, with most of the correction due to the diagrams shown in 
figure 2; this spoils the perturbation expansion. Note that this is not a problem 
when rriH- > rrit, as in that case the kinematics of the final state H~b system forbids 
the intermediate i from going on-shell. 

However, as is clear from the figure and the fact that the problem occurs when 
the i is on-shell, the diagrams can also be thought of as the production of a top 
quark pair at leading order in QCD, with decay of the t particle producing a charged 
Higgs boson. That is, top pair and single top production interfere beyond LO in the 
single top cross-section, and the question arises of whether it is possible to separate 
the two processes, and thus maintain a meaningful definition of H~t production at 
higher orders in perturbation theory. It must be stressed that such a separation 
is an approximation for practical purposes only (i.e. H~t and tt production really 
do interfere at the quantum level), and is dependent on the experimental analysis 
cuts which are applied. However, and as we will see, such a separation is indeed 
possible subject to cuts, allowing the single and pair production processes to be 
added incoherently to a sufficiently good approximation. 

There are a number of advantages that result from separating H~t from top 
pair production, in trying to accurately represent their sum. A theoretically rigorous 
superposition of single and top pair production would require at that LO top 

pair production (with decay of the i) be added to NLO H~t at the amplitude level. 
Such a procedure includes all interference effects, but neglects NLO contributions to 
the top pair process. The latter are known only in the narrow width approximation 
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in which the t quark is produced on-shell, and are known to be large (of order 50% for 
cuts used to isolate the top pair production cross-section). If one instead combines 
the single and top pair processes incoherently, the NLO corrections are included at 
the expense of the interference term. Importantly, for signal cuts used to isolate the 
single top process, the NLO corrections to top pair production are, in a well-defined 
sense, larger than the interference between the single and pair production processes. 
The most accurate description then results from adding the H't and top pair matrix 
elements incoherently, thereby neglecting the interference term. 

Such a discussion is not meaningful unless one has a means of quantifying the size 
of the interference between single top and top pair production. If this is indeed the 

case, then the systematic uncertainty due to interference effects can be compared 
with other uncertainties in a given analysis (such as scale variation), in order to 
ascertain whether or not it is legitimate to regard H~t and ti — > tH~b as incoherent 
processes. 

The interference problem described above is by no means restricted to H~t pro- 
duction, but also occurs in other contexts such as single top quark production in 
association with a W boson. Solutions to the problem in that context have been 
widely studied [33, 34, 35, 36]. In [10] two definitions of the Wt process were given 
which could be applied in a fully exclusive context i.e. within a parton shower frame- 
work. These definitions were called Diagram Removal (DR) and Diagram Subtrac- 
tion (DS), and were subsequently implemented in MC@NLO. The definitions were 
designed so that the difference between them measured the degree of interference 
between Wt and tt production, where the t in the latter case decays to a Wb pair. 
Thus, the difference between the DR and DS results can be used to estimate the 
systematic uncertainty due to interference effects. This will not be small in general, 
and depends on the analysis cuts applied to the final state. 

A phenomenological analysis was carried out using these results in [32], which 
showed that for signal cuts used to isolate the Wt signal, the interference with top pair 
production is indeed small with respect to other uncertainties in the analysis (most 
notably, the scale variation of the cross-section which is representative of the size of 
higher order corrections). This was also the case when Wt and top pair production 
were themselves backgrounds to a third process (the specific example of Standard 
Model Higgs production followed by decay toaW boson pair was considered) . In such 
cases, the most accurate description of the relevant single and top pair production 
processes results from an incoherent sum of the two, for which one can use MC@NLO 
as detailed in those papers. The strong theoretical similarities between the Wt and 
H^t amplitudes imply that the same definitions can be applied in the latter case, 
and that they will be similarly successful from the phenomenological point of view. 
This is indeed the case, as we will see in the following subsections. First, we recap 
the definitions of DR and DS, with the discussion tailored to the present context. 
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3.1 Diagram Removal (DR) and Diagram Subtraction (DS) 

In DR, the diagrams of figure 2 (which we may call doubly resonant due to the 
intermediate top quark pair) are removed at the amplitude level. Thus, the squared 
NLO amplitude for H~t production has no interference with top pair production by 
construction. A potential difficulty of this approach is that both electrowcak and 
QCD gauge invariance are violated. However, this was seen not to be a problem in 
practice (see [10] for a detailed discussion). 

In DS, the cross-section for H~t production is modified by a local subtraction 
term, which effectively removes the contribution of doubly resonant diagrams point- 
by-point in phase space. This is a gauge-invariant procedure, although there is a 
degree of arbitrariness in the explicit construction of the subtraction term. Schemat- 
ically, one has 

da^-f = daH-t — da^-^, (3.1) 

where the first term on the right-hand-side denotes the fully exclusive cross-section 
for H~t production (including doubly resonant diagrams), and the second the sub- 
traction term. The latter must satisfy the following requirements: 

1. When the invariant mass of the final state H~b system satisfies ruu-i = mt, 
the subtraction term must be equal to the fully exclusive cross-section for top 
pair production, with the t decaying to H~b. 

2. The subtraction term must fall away sharply as mjj-f, moves away from the 
top mass. 

The first condition amounts to subtracting the doubly resonant contribution as re- 
quired, and the second ensures that the diagrams in figure 2 contribute unmodified 
when the t is off-shell. These are the only two conditions that the subtraction term 
must satisfy, and following [10] we use the explicit form 



(3.2) 



fswimt) 

Here fsw is the Breit-Wigner function, and A^*^ denotes the LO scattering amphtude 
for tH~b production including only doubly resonant diagrams (note that these form 
a gauge-invariant set), but with the kinematics reshuffled so as to place the t on- 
shell. We obtain the relevant matrix elements from MadGraph [37, 38]. The use 
of the full amplitude ensures that spin correlations of the i arc present in the decay 
products H^b (these arc not to be confused with the decay products of the top quark 
and charged Higgs boson, which do not currently have spin correlations included). 
Clearly this squared amplitude fulfills the first of the above requirements. It is then 
multiplied by a ratio of Breit-Wigner functions, which damp the subtraction term 
as moves away from the top mass, thus fulfilling the second condition above. 

More details can be found in [10]. 
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Given the fact that the subtraction takes place at the cross-section level in DS, 
rather than at the amplitude level as in DR, the interference term between H^t and 
tt is still present in DS. It follows that the difference between DR and DS mostly 
measures the interference, as stated above. The difference between DR and DS is also 
affected by ambiguities in the formulation of the subtraction term, and by potential 
gauge invariance violation in DR. However, both of these effects are small [10]. 

We will see in section 4.2 that the difference between DR and DS can indeed be 
made small subject to adequate cuts. Firstly, however, we present example results 
for the large Iliggs mass region mn- > irit, where the complications of this section 
are not relevant. 

4. Results 

In this section, we present example results from our MC@NLO implementation for 
H^t production. Our aim is not to undergo a thorough phenomenological analysis, 
but rather to present a few observables which demonstrate the differences between a 
fixed order description, and one that is matched to a parton shower. In the following 
results, we consider the LHC at 14 TeV, set renormalization and factorization scales 
to IJ>F — IJ'R — {ttih- + "^t)/2, and use a top quark mass of rrit — 172 GeV. For 
the charged Higgs Yukawa couplings we assume a type-II two-Higgs doublet model 
(eq. (A. 14)), with tan/5 = 30 and the NLO running bottom and top quark masses 
with scales as above. However, tan f3 is simply a scaling factor of the bg tH~ cross 
section and once it lies above tan /9 > 10 is has no impact on the relative size of the 
NLO corrections or the distributions of the final state particles [13, 14, 21]. We use 
CTEQ5M1 partons throughout [39]. Although more recent global parton analyses are 
available (e.g. [40, 41]), this choice facilitates comparison with the results of [14, 21] 
where appropriate. 

4.1 Lcirge chcirged Higgs mass: niH- > rrit 

In this section we consider, unless otherwise stated, mH-—^00 GeV. We begin by 
showing the transverse momentum and rapidity distributions of the top quark and 
charged Higgs boson in figure 3. One sees that there is not a great deal of difference 
between NLO and MC@NLO descriptions of these observables, as perhaps expected 
given the inclusive nature of these observables, and the fact that they are not sub- 
ject to logarithmic corrections. Instead, the pt and rapidity distributions act as a 
consistency check between the NLO and MCCilNLO calculations. 

We now consider observables designed to manifest the differences between the 
fixed order and parton showered approaches. The quantity 

Pt,sum = \Pt,H + Pt,top\ (4.1) 
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Figure 3: Comparison of NLO and MC@NLO results, with parameters as given in the 
text. Shown are the transverse momentum and rapidity distributions of the top quark 
(upper hne) and Higgs boson (middle line), the transverse momentum of the H^t system 
(bottom left) and the azimuthal angle between the top quark and charged Higgs boson 
(bottom right). 



(i.e. the magnitude of the vector sum of transverse momentum of the H~t system) 
is shown in figure 3 (bottom left). At LO, this distribution would be a Dirac delta 
function at Pt,sum = due to 4-momentum conservation. At NLO, there is a sharp 
rise as pt,sum due to the real emission contributions, whilst the zero bin is 
negative due to virtual corrections. This behavior is smoothed by the parton shower, 
and the MC@NLO curve indeed displays the characteristic Sudakov behavior. 
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Also of interest is the distribution of the azimuthal angle between the top quark 
and charged Higgs, shown in figure 3 (bottom right). At LO, this would be a delta 
function at = tt i.e. the final state particles are produced back-to-back. At NLO, 
this gets decorrelated by the emission of one hard parton. There is also a suppression 
of the (f) — TT bin due to the virtual corrections, resulting in a slightly peculiar shape in 
the last two bins. When the full parton shower is added this feature is smoothed out 
to produce a more physical-looking decorrelation, thus demonstrating the advantages 
of the parton shower approach. 

The results of this section indicate that the MC@)NLO is indeed working as 
expected. We now consider some phenomenological properties of b and light jets. Our 
motivation is as follows. Previous analyses of the H~t process (such as [42, 43, 44]) 
have argued that additional b jets, by which we mean b jets that are not the hardest 
b jet and thus are less likely to have come from the decay of the top quark, have 
different kinematic properties (specifically the transverse momentum distribution) to 
light jet radiation. Furthermore, that one may exploit this difference to design event 
selection criteria, i.e. reduce backgrounds. A typical example is the background 
due to vector boson plus multijet (or pure multijet) production. In a non-negligible 
fraction of cases, light jets in the multijet background may be mistagged as b jets, 
and therefore be mistaken for signal events. It is much less likely, however, that two 
light jets will be mistagged, such that if one can convincingly identify events with 
two b jets, the signal to background ratio may be significantly enhanced. 

This naturally raises two questions. First of all, what are the fractions of events 
in which radiated light jets or additional b jets are present? Secondly, are additional 
b jets harder than radiated hght jets? An MC@NLO description is better suited to 
answering this question than a purely fixed order description, owing to the greater 
number of final state partons present, and we therefore consider these questions here. 

In the following we will consider cases where the top quark decays both leptoni- 
cally, and also hadronically. The former case is more promising experimentally, due 
to the absence of hard jets from the top decay (we will see this in more detail in 
what follows). We will assume, unless otherwise stated, 100% 6-tagging efficiency. 
Furthermore, we will not consider systematic uncertainties due to reconstruction of 
the top quark or charged Higgs kinematics (i.e. we will assume that the decay prod- 
ucts are sufficiently hard in pt to be observable). We cluster jets according to a kt 
algorithm with D — 0.7, and require all jets to lie in the 6-tagging detector volume 

\ri\ < 2.5 Pt > 25 GeV. (4.2) 

In any given event, one sees a number (possibly zero) of 6 jets, but has no way of 
knowing which b jet came from the decay of the top quark. Thus, in comparing the 
properties of additional b jets with radiated light jets, one assumes that the hardest 
b jet in any event arises from the top quark, and that the second hardest b jet arises 
from QCD radiation (in our calculation this happens either in the NLO matrix 
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element, or in the parton shower; hkewise for hght jets that do not come from a top 
decay). Figures 4 and 5 show the transverse momentum and rapidity distributions 
of the two hardest bottom and hght-flavor jets, if present, for two charged Higgs 
masses 300 GeV and 800 GeV. The hght jet results are shown for leptonic as well 
as hadronic top decays. In the first row of distributions we see the behavior of the 
two bottom jets, one from the top decay and the other (mostly) from initial-state 
gluon splitting. The harder of the two bottom jets peaks at transverse momenta 
around 50 GeV, as indeed is expected from the three-body decay of a top quark. 
The position of this peak therefore does not change with the Higgs mass. The softer 
of the two bottom jets is peaked at small transverse momenta and extends to larger 
rapidities than the decay jet — the typical pattern of initial state radiation. Thus, 
the assumption that the first hardest jet arises from the top decay and the second 
hardest from QCD radiation appears to be well-founded. At the higher Higgs mass 
(i.e. a larger scale in the hard process) both rapidity distributions become flatter. 
This can happen because of enhanced coUinear radiation in the presence of a harder 
QCD scale, or because of momentum conservation effects between the harder Higgs 
boson and its accompanying jets. 

For leptonic top decays all light-flavor jets arise from QCD radiation, rather than 
from decay of the top quark. This is reflected in the sharp drop of the distributions 
as well as the flat r] behavior. As mentioned above, it is interesting to compare the 
properties of the second hardest b jet (i.e. the b jet mostly coming from QCD 
radiation rather than the top quark decay) with the hardest light jet. In particular, 
one may compare the transverse momentum distributions shown in figure 4. One sees 
that the difference between the second hardest b jet and hardest light jet distributions 
is not large. If anything, for ttih- = 300 GeV the pr distribution of the hardest light- 
flavor jet tends to fall off slower than the second hardest bottom jet distribution. This 
behavior is more clearly visible for heavy Higgs masses, where the px distribution 
reaches half its maximum value at 55 GeV for the second bottom jets and at 65 GeV 
for the hardest light-flavor jet. In principle, both the light-flavor jet and the bottom 
jet are radiated coUinearly, so the source of this slight discrepancy is simply that 
there are more light-flavor jets to pick from. 

One may also consider hadronic decays, and results for the transverse momen- 
tum and rapidity distributions of the two hardest light jets are shown in the lower 
lines of figures 4 and 5. One sees that the hardest jets become significantly harder 
and more central, as expected from the fact that they now come predominantly from 
the top decay rather than QCD radiation. Unfortunately, this distinguishing feature 
is limited in its experimental use, because the high mass scales in this process (stem- 
ming from the top quark and Higgs boson) allow for initial and final state radiation 
comparable in hardness to the jets from the top decay [45, 46, 47]. Indeed W bosons 
and top quarks that undergo hadronic decays can only be reconstructed in special 
kinematic regions, such as when they are highly boosted [48, 49, 50, 51, 52]. 



16 





Figure 4: Normalized transverse momentum distributions of the two hardest jets. We 
include the detector cuts of eq. (4.2). First row: the two hardest b jets; second row: the 
two hardest light-flavor jets for a leptonic top decay; third row: the two hardest light-flavor 
jets for a hadronic top decay. The left-hand column corresponds to rriff- = 300 GeV, and 
the right-hand column to uih- = 800 GeV. 

Above, we have compared whether additional b jets are harder than radiated 
light jets. One must also consider (in light of the first question above) the likelihood 
of observing events with radiated light jets or more than one b jet, subject to the 
detector constraints. 

In panels (a) of Table 1 we show the probability to see a light-flavor jet in addition 
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Figure 5: Normalized rapidity distributions of the two hardest jets. We include the 
detector cuts of eq. (4.2). First row: the two hardest 6 jets; second row: the two hardest 
light-flavor jets for a leptonic top decay; third row: the two hardest light-flavor jets for 
a hadronic top decay. The left-hand column corresponds to mu- = 300 GeV, and the 
right-hand column to m^/- = 800 GeV. 

to the bottom jet from the top decay, for a number of different cases of the detector 
cuts 

\v\< Vcut Pt > Pt,cmu (4.3) 

and for leptonic decays of the top quark (i.e. where light jets can only come from QCD 
radiation). As expected, such a light QCD jet from soft and coUinear jet radiation 
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45 GeV 


83.2 


79.2 


72.3 


61.0 


39.9 


87.1 


83.3 


76.8 


65.7 


45.2 


65 GeV 


60.9 


57.3 


51.7 


43.2 


28.8 


70.5 


66.9 


61.3 


51.9 


35.9 




85 GcV 


44.4 


41.5 


37.1 


31.1 


21.3 


5(3.2 


53.3 


48.0 


41.0 


28.7 




25 GeV 


17.8 


14.3 


10.0 


5.7 


2.3 


21.4 


17.1 


12.1 


7.1 


3.0 


(c) 


45 GeV 


12.9 


10.6 


7.6 


4.5 


1.8 


16.6 


13.7 


9.9 


6.0 


2.5 


65 GeV 


9.4 


8.0 


5.9 


3.5 


1.6 


12.3 


10.5 


7.9 


4.9 


2.0 




85 GeV 


7.2 


6.4 


4.8 


3.0 


1.4 


9.7 


8.5 


6.5 


4.1 


1.7 



Table 1: Probability (%) to see (a) a light jet given that a b jet has been observed, where 
the top decays leptonically; (b) a light jet given that a b jet has been observed, where the 
top decays hadronically; (c) an additional b jet given that one b jet has been observed, 
where the top decays leptonically. All jets satisfy eq. (4.3). 

appears in a sizeable fraction of the events. The lower percentages compared to 
heavy sgluon or squark/gluino production [45, 46, 47] are due to the fact that here 
we only consider central jets out to \ri\ = 2.5. This percentage increases with the 
hard scale in the process, such that for a heavy Higgs more than half of all events 
show a central jet within I77I < 2.5 and above 25 GeV. Even above the typical range 
of jets from W decays, i.e. requiring px > 45 GeV, the fraction of events with at 
least one QCD jet varies between 32% and 42%, depending on the Higgs mass. 

In panel (b) of Table 1 we show similar results, but for hadronic top decays, 
where light jets may arise either from QCD radiation or from the decay of the W 
boson from the top quark. This increases the average number of light-flavor jets 
per event by (roughly) two. Indeed, in panel (b) of Table 1 we indeed see that the 
percentage of events with jets above 25 GeV reaches close to 100%. 

In panel (c) of Table 1 we show the probability of observing an additional b 
jet, given that a hard b jet has already observed. Again, as expected from collinear 
radiation the percentage of events with two bottom jets decreases sharply with more 
central or harder bottom jets. Only in the most central phase space region with 
Pt ~ 25 GeV does the probability of observing two bottom jets reach 0(20%). 
The shghtly higher probabilities for larger Higgs masses is an effect of the generally 
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hardened Pt spectrum of the jet radiation. Considering the full detector region of 
eq. (4.2), the probability of seeing a second b jet given a first is 17.8% (21.4%) for 
a light (heavy) charged Higgs boson. Again, from a QCD perspective this overall 
percentage is very low [45, 46, 47] because we require \r)\ < 2.5 to allow for b tagging, 
instead of the full jet range \ri\ < 4.5. However, this prediction is overly optimistic in 
the sense that we have assumed 100% 6-tagging efficiency. The numbers in Table 1 
still have to be multiphed with a b tagging probabihty around 50% — 60% (or indeed 
even lower for the smaller pt,cut values), which means that CMS and ATLAS can 
expect of the order of 10% of their H~t events to include a second tagged bottom 
jet, if a first has been seen. This number is subject to higher order corrections to the 
gg — > btH^ process [53, 54], but these do not change things significantly. 

The upshot of the above discussion is that one can observe two b jets in only 
a small fraction of events. Furthermore, comparison of the transverse momentum 
distributions seems to indicate that additional b jets are not substantially harder 
than radiated light jets. One may make this latter point more specific as follows. 
Consider that a hard b jet has been observed. We may then ask the question: 
what is the probability that the hardest jet is a 6 jet? This is shown in panels (a) 
and (b) of Table 2 for the cases where the top decays leptonically and hadronically 
respectively. For leptonic decays, the probability is reasonably high (i.e. upwards 
of 80%), reflecting the fact that for leptonic decays there is a hard b jet from the 
top decay, but light jets arise only from QCD radiation. However, it is interesting 
(and perhaps surprising) to note that for pT,cut ~ 25 GeV and 77cut — 2.5, a radiated 
light jet is harder than the primary b jet in a sizeable fraction of cases (i.e. around 
20%). For hadronic decays (panel (b)), the probability that the hardest jet is a 6 jet 
is much lower, given that there is competition from the light jets from the top decay. 
The numbers change only slightly in going to a higher Higgs boson mass. 

Given that, as remarked above, one really wants to isolate two b jets in order 
to reduce multijet backgrounds, the more experimentally relevant question is: given 
that one hard b jet has been observed, how likely is it that we find the second b 
jet by asking for the two hardest jets in the event? This probability is shown in 
panels (c) and (d) of Table 2, and is considerably below 50% even for the best-suited 
phase space regions. Thus, for central jets in the case of leptonic top decays we 
see that the radiated light-flavor jets are typically harder than the additional b jet. 
This is simply a combinatorial effect — from collinear radiation we expect either one 
or three bottom jets, where three bottom jets are strongly suppressed. For regular 
QCD radiation we have many more diagrams leading to additional jets, so when 
we ask for the hardest radiated jet it is usually one bottom jet vs. the hardest of 
several light-flavor jets. As expected, this fraction drops sharply for the (already 
experimentally less promising) hadronic top decay signature, as can be seen in panel 
(d) of the table. 

The above results demonstrate that isolating events in which there are two hard b 
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9.0 


7.1 


14.5 


13.8 


12.1 


9.9 


7.2 


65 GeV 


17.1 


16.2 


14.2 


10.9 


7.8 


16.8 


16.0 


13.9 


10.9 


7.5 




85 GeV 


20.1 


19.3 


17.3 


13.3 


8.9 


18.4 


17.6 


15.3 


11.9 


7.8 



Table 2: Probability (%) that (a) the hardest jet is a 6 jet, given that at least one b jet 
has been observed, where the top decays leptonically; (b) the hardest jet is a 6 jet, given 
that at least one b jet has been observed, where the top decays hadronically; (c) probability 
that the two hardest jets are both b jets, where the top decays leptonically; (d) probability 
that the two hardest jets are both b jets, where the top decays hadronically. 

jets is difficult. There are not many events in which two or more observable b jets are 
present, and their transverse momentum properties do not significantly distinguish 
additional b jets from radiated light jets, even in the case of leptonic decays (when 
there are no hard light jets from the top quark decay). The above results, however, do 
not constitute a thorough phenomenological analysis, such that further investigation 
may be useful. 

The preceding discussion is restricted to the region of large Higgs mass, where 
one does not have to worry about interference efi^ects. One must also consider the 
regime of small Higgs mass, which is the subject of the next section. We will focus 
on interference issues, rather than properties of b and light jets. However, it is worth 
bearing in mind that the above discussion would be significantly modified for small 
Higgs mass, where the second 6 jet may indeed prove to be useful in reducing top 
pair production backgrounds (analogously to the Wt process studied in [32]). 
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4.2 Small chcirged Higgs rriciss: niH- < rrit 



As described in section 3, when the charged Higgs boson mass is lower than the 
top mass, the H~t process at NLO interferes with top pair production at LO, with 
decay of the t particle into a charged Higgs boson and (anti)-bottom quark. It is still 
possible to construct an MC@NLO for H^t production, using the diagram removal 
and diagram subtraction definitions, first applied in the context of Wt production 
in [10]. MC@NLO then gives a sensible result only subject to the approximation that 
the single and top pair processes can be added incoherently i.e. that the interference 
between the two processes can be neglected. This will not be true in general, but 
is likely to be true for analysis cuts used to increase the ratio of the single and top 
pair cross-sections. The DR and DS MC@NLO codes are defined independently of 
any subsequent analysis cuts, and the difference between these codes can be used to 
estimate the systematic uncertainty due to interference effects. 

The object of this section is to both clarify and strengthen the above remarks, 
by comparing the DR and DS results for various charged Higgs boson masses. In 
order to verify that the DR and DS codes function as required, it is useful to have 
a continuous cut parameter that smoothly interpolates between the regions where 
the interference is expected to be small and large respectively. Following [10], we 
introduce a transverse momentum veto on the second hardest B hadron occurring in 
an event. That is, one first finds all B hadrons satisfying the detector constraint 



where rj is the pseudo-rapidity. Then events are rejected if a second hardest B hadron 
is present, with transverse momentum satisfying: 



The reasoning is that top pair-like events at the LO parton level contain two b quarks 
in the final state, whereas H~t-\ike events contain only one. At NLO plus parton 
shower level, one still expects top pair-like events to contain roughly two observable 
B hadrons, and single top-like events to have less than two. Thus, for small values 
of the Pt,veto; the top pair process is efficiently reduced relative to the single top 
mode. Large values of the veto (pt.vcto oo) correspond to no constraint, and an 
interference that is not expected to be small in general. 

Note that there are different choices in how to apply this veto. One can apply it 
even in the pure NLO calculation, by using b quarks rather than B hadrons, although 
there are difficulties interpreting this phenomenologically (see the detailed discussion 
in [10]). In the parton shower approach, one could also apply the veto at the level of 
b jets for example. However, the aim here is merely to check the mutual validity of 
the DR and DS procedures, rather than to be phenomenologically complete. Also, 



V\ < 2.5, 



(4.4) 



Pt > Pt.veto- 



(4.5) 
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Table 3: Total MC@NLO cross-section results (in pb) for H~t cross-sections subject to 
the transverse momentum veto of eq. (4.5), with parameters as described in the text. Both 
the DR and DS results are shown, and errors shown are due to variation of the common 
renormalization and factorization scale by a factor of two. The statistical error on the ratio 
values is less than a percent. 

applying the veto in the same way as in the Wt case of [10] allows one to compare 
the size of interference effects in both cases. 

For renormalization / factorization scale choices and top quark mass fixed as 
above, the level of disagreement between DR and DS will depend in general on the 
charged Higgs boson mass as well as the transverse momentum veto. Thus, we 
evaluate cross-sections for a number of choices of mu- < mt- Results are shown in 
table 3, also for a number of choices of pt,veto- The quoted uncertainties are due to 
variation of the common renormalization and factorization scale by a factor of two. 
This uncertainty can be larger if the scales are allowed to vary independently. One 
sees that the level of agreement between DR and DS is around 4-5% at Pt,veto — 10 
GeV. This is higher than in the case of Wt production, indicating that interference 
effects are larger in the present case (they also depend upon the Higgs mass of course). 
Nevertheless, the DR and DS results overlap within the scale variation uncertainty. 
The results with no veto applied show a marked disagreement (greater than 10% in 
some cases). This shows that interference with tt production is indeed a problem in 
some phase space regions, as is to be expected. Interestingly, for the value of the 
charged Higgs mass that hes closest to the top quark mass, the results agree well 
even when no veto is applied. We will discuss the threshold region in more detail in 
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the following sectfon. 

Total cross-sections are only part of the story. One must also check that inter- 
ference effects are small locally in phase space, and this is possible by comparing 
kinematic distributions obtained from both the DR and DS codes. Examples are 
shown in figure 6 for pt,veto = 10 GeV, and mn- — 100 GeV (results are repre- 
sentative of other charged Higgs masses). One indeed sees good agreement in the 
observables considered. Some comments are in order. Firstly, we have here shown 
only a few observables, with a particular cut designed to reduce the interference with 
top pair production (i.e. the transverse momentum veto). A more realistic analysis 
designed to isolate the charged Higgs production process would use different cuts 
to reduce the interference, such as jet vetoes (sec e.g. the Wt analysis of [32]). It 
may still be possible to define observables in such cases that are sensitive to inter- 
ference effects, and whether or not this is the case can be ascertained by comparing 
the output of the DS and DR predictions. Nevertheless, the results here indicate 
that the incoherence of single and top pair-instigated charged Higgs production is 
a reasonable approximation, over a sufficiently large region of the phase space for 
MC@NLO to provide an accurate description when possible (as was found to be the 
case in Wt production in [32]). In the next section, we explain how to interpret the 
behavior of the MC@NLO calculation across the threshold mn- = rrit. 

4.3 Threshold behavior 

In the previous sections, we have outlined the calculation of the H~t production pro- 
cess up to NLO, and described the interface to a parton shower with the MC@NLO 
formalism. For mn- > TUt this is straightforward, whereas for mu- < rrit a problem 
arises due to interference with top pair production. The aim of this section is to 
clarify any possible confusion in the behavior of the cross-section across the top mass 
threshold. 

A consequence of the subtraction inherent in the DR and DS definitions is that 
the MC@NLO prediction for H~t itself is not continuous across the threshold ttih- — 
rrit. Nor should it be, given that the sum of the single and top pair processes is 
continuous, and the top pair contribution is subtracted out in either DR or DS. The 
result is that there is a discontinuity in the H~t calculation (defined using either 
DR or DS) as ttih- tends towards rrit from below. This discontinuity would (in 
MC@NLO) be filled in by the addition of top pair contributions, which in any given 
analysis must be added background. 

This is depicted explicitly in figure 7, which shows the fully inclusive total H^t 
cross-section calculated using DR and DS as a function of the charged Higgs boson 
mass, with all other parameters set as above (in particular rrit = 172 GeV). The 
curves disagree for the low mass region, as expected if no analysis cuts are applied to 
reduce the interference with top pair production. The apparently unphysical behavior 
near rrin- — rrit is purely a consequence of the fact that MC@NLO is calculating a 
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Figure 6: Comparison of DR and DS results, for mjj- = 100 GeV, and other parameters 

as outlined in the text. Shown arc the transverse momentum and rapidity distributions of 
the top quark (upper line) and Higgs boson (second line), as well as the azimuthal angle 
between the top quark and charged Higgs boson (lower line). Error bars denote statistical 
uncertainties. 

qualitatively different cross-section on either side of this value. For the high Higgs 
mass region ttih- > mt, the curves agree due to the fact that no subtraction is 
apphed. 

We stress that no interpolation is necessary, or indeed meaningful, between the 
calculations on either side of the mass threshold. Furthermore, the discontinuity is 
not problematic in any real analysis, in which the H~t prediction must be combined 
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Figure 7: The total NLO cross-section for H~t production, shown for DR (sohd) and DS 
(dashed), for a top mass of 172 GeV. The error bars correspond to statistical uncertainties. 
The discontinuity is due to subtraction of resonant contributions for mjj- < mt, and would 
be filled in after adding in the top pair production background. 



with the corresponding top pair process (a discussion of threshold behavior can also 
be found in the context of a strictly NLO calculation in [21], where however the DR 
and DS definitions have not been used). 

5. Conclusion 

In this paper, we have considered the process of charged Higgs boson production in 
association with a top quark. This is an important scattering process at the LHC, 
given that charged Higgs bosons generically occur in extensions to the Standard 
Model. In order to increase the accuracy with which this process can be calculated, 
we have implemented it in the MC@NLO framework for combining fixed order matrix 
elements at next-to-leading order with a parton shower. 

The details of the MC@NLO calculation are different in the two kinematic regions 
in which the charged Higgs boson mass is less than and greater than the top mass 
respectively. For mn- > mt, the NLO calculation of the H~t process is well-defined, 
and the implementation in MC@NLO follows the procedure adopted in other single 
top (and indeed non-top related) processes [8, 26, 9]. 

In this high Higgs mass region, wc presented example kinematic distributions 
involving the top and charged Higgs bosons, which demonstrated the expected dif- 
ferences between the fixed order and parton showered approaches. We also looked 
at the properties of b jets in addition to the top decay, to see if these are in any way 
different to radiated light jets (motivated by previous analyses which aim to exploit 
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such differences). The conclusion is that radiated b jets are not substantially harder 
than light jets. Furthermore, the fraction of events containing two b jets is small. 

For the low Higgs mass region niH- < rrit, the NLO cross-section for H^t pro- 
duction becomes ill-defined due to interference with top pair production at LO, with 
decay of the t to a charged Higgs boson and b quark. There is then some ambi- 
guity in how to proceed, and we proceed by analogy with the Wt case discussed 
in [10, 32]. That is, one may approximate the sum of the single and top pair produc- 
tion processes by an incoherent sum, subject to adequate cuts on the phase space. 
We produced two MC@NLO implementations using the diagram subtraction (DS) 
and diagram removal (DR) approaches, such that the difference between the two 
codes provides an estimate of the systematic uncertainty due to interference effects. 
The resulting software can then be used in the low Higgs mass region to efficiently 
generate H~t-like events. 

Using the DR and DS codes we find that, whilst the interference is sizeable 
when no cuts are applied (as is to be expected), the DR and DS definitions agree 
well when a suitable cut is applied. Here we used the simple example of a transverse 
momentum veto on the second hardest B hadron in an event. The results are similar 
to the Wt case of [10] (albeit with a larger interference effect seen in the present 
case), and therefore suggest that isolation of the charged Higgs signal is possible 
in more realistic analyses (based on e.g. b jet vetoes). We postpone more detailed 
phenomenology to future work. 
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A. Renormalization of Yukawa couplings 

In this appendix we derive the counterterms resulting from the renormalization of 
the charged Higgs Yukawa couplings. We first examine the case of a scalar neutral 
Higgs boson, before considering the example of a type II two-Higgs-doublet model. 
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Similar remarks would apply in the case of a type I model. Our presentation follows 
that of the neutral Higgs boson case presented in [55]. 

The Yukawa coupling of the Higgs boson to quarks in the Standard Model is 
given by: 

>CYuk. = -^K^ + H)QlQr, (A.l) 

where v is the vacuum expectation value of the Higgs boson field H, and Ql,r are 
the left and right-handed quark fields. Prom the above Lagrangian one sees that the 
bare quark mass and Yukawa couplings are related by: 

m. = ^. (A.2) 

At one-loop order the bare mass and Yukawa couplings are renormalized according 
to: 

Ao = /x^A(l + 5A); (A.3) 
mo — m + Sm, (A. 4) 

where 5X and Sm are the appropriate counterterms, and A and m (scheme-dependent) 
renormalized quantities. One may choose a renormalization scheme such that the 
Higgs VEV V is not renormalized at one loop order. Then one finds 

SX = — , A.5 

m 

i.e. that the Yukawa coupling renormalization is simply related to the mass renor- 
malization. Note that 5m oc m, such that the Yukawa coupling counterterm is 
independent of the quark mass. In the MS scheme, this has the form 

■ MS Cpas f 3 



^^Ms ^ \ ~ 31n47rj . (A.6) 

We now consider the case when two Higgs doublets are present, and focus on the 
case of a type II two-Higgs-doublet model. Focussing on one quark generation, the 
coupling of the Higgs doublets to the fermions is given by (see e.g. [6]) 

^yT" = -^ij [(Aft + 5X,)hnH{Q{ + (A, + 5\)tRQ{Hi + ...], (A.7) 

where the ellipsis denotes terms obtained from those shown by Hermitian conjuga- 
tion. Here R and L denote right- and left-handed quarks respectively, and i and j 
are flavor indices. The Higgs and fermion SU(2) doublets are given by 

H^=[fX ^^=(^0 1' Q^=l r h (A-8) 
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where and are the neutral and charged Higgs boson respectively. Expanding 
each Higgs doublet about its vacuum expectation value corresponds to the replace- 
ments (including conventional factors of \/2) 

where P° correspond to the pseudo-scalar Higgs particles and Goldstone bosons. 
Substituting eq. (A.9) into eq. (A. 7), one may then apply an analogous argument to 
the Standard Model case to each Higgs doublet separately, to obtain the relations 

= a, = ^. (A.10) 

It is conventional to define the parameters v and fi parameter from the two Higgs 
VEVs i.e. 

^ f 2 f sin /3 / . . , N 

tan/3 = — = ^. A.ll 

vi V cos p 

Then the coupling of the physical charged Higgs field to fermions can be written 

rrib tan/3(l - 75) + rrit cot (3{1 + 75)] , (A.12) 



where the appropriate CKM matrix element has been factored out. This has the 
form 

GH-tb = ^Vtb{a-lrf,), (A.13) 

where 

a = (m;,tan/3 -|- cot/3) , 6 = — p (m^ tan /3 — cot /3) . (^-14) 
V 2 v 2 

Prom the above remarks, these clearly have the corresponding one-loop renormaliza- 
tion 

ao^a(l + da), bo = b(l + Sb), 5a = db = — , (A.15) 

m 

where m is the mass of either the top or down-type quark, given the independence 
of the final result on the mass. 



B. Calculation of M 

The counterterms to real-emission matrix elements in subtraction formalisms have 
to cancel locally the relevant divergences, in order for the numerical integration to be 
well-defined and stable. The case of the coUinear branching of a gluon (in the initial 
or final state, with the gluon being spacelike and timelike respectively) merits some 
attention because it has a non-trivial local azimuthal dependence, that vanishes when 
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the azimuthal integration is carried out. As shown in ref. [18], in the FKS subtraction 
formahsm the azimuthal-dependent part of the local collinear counterterm is quite 
simple, and is formally identical to the better-known azimuthal-independent one: 
both have the structure of a kernel (that depends only on the parton identities, but 
not on the hard process), times a short-distance. Born- like cross section. 
To be more specific, consider the process 



+ P{P2) X{K) + S{k), (B.l) 

where Greek letters represent (massless) partons, and X any other particles that 
may be present (i.e. in the present case, X = t + H~). In the collinear limit p2 || 
the squared matrix element for the above process is given by (see eq. (B.41) of [18]) 

M(pi,P2) ^ ^ \p{z)M^'^\p,,zp2)+Q{z)M{puzp2)\ , (B.2) 
k-p2 L J 

where P are the usual Altarelli-Parisi kernels, Q are other universal kernels (given 
at the leading order in eqs. (B.42)-(B.45) of ref. [18] for initial-state collinear split- 
tings, and in eqs. (B.31)-(B.34) of that paper for final-state collinear splittings; these 
kernels arc thus different for spacclikc and timclikc branchings already at the lead- 
ing order), is the relevant Born contribution, and M. is a, Born- like function, 
which however keeps track of the azimuthal correlations in the branching process. 
The contribution of QM vanishes upon integration over the azimuthal angle of the 
branching, which is why this term can be neglected in the analytical computation of 
the collinear divergences in 4 — 2e dimensions. Locally, it is different from zero if the 
parton involved in the branching which enters the hard reaction is a gluon, and there- 
fore needs to be taken into account for the construction of the local counterterms in 
a (efficient) numerical NLO program. 

The helicity interference amplitude is given by 

^ = ^Re|gH|^fU(o)|^ 



and in the case of Wt production in [10], we calculated this in the helicity formalism. 
Here we adopt a different (and somewhat quicker) approach, involving projection of 
the helicity dependent Born amplitude with relevant tensors. 

Let M^^^ be the Born amplitude of figure 1 before contraction with the gluon 
polarization tensor {—g^^)- Then the helicity interference is given by [56] 
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where k± is the transverse component of the gluon momentum entering eq. (B.l). 
Carrying out the calculation in the present case gives 



4Ar^(7V2 _ 1) 

4(2cos^0- l)(m( — rn\j-){ml{m\^ - u) -\- u{—m'\j- -\- s -\- u)) 

s{m1 — u)' 



^ 17112 7x2 ' y^-"^) 



where (f) is the gluon azimuthal angle, and we have defined the Mandelstam invariants 
s = (pi+P2)^, ti^t—ml^{ki—pif—ml, Ui^u—ml^{k2—pif—ml. (B.6) 
One sees immediately that this term vanishes upon integration over 0. 
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